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Abstract 

We numerically integrate the semiclassical equations of motion for spher- 
ically symmetric Einstein-Maxwell theory with a dilaton coupled scalar field 
and look for zero temperature configurations. The solution we find is studied 
in detail close to the horizon and comparison is made with the corresponding 
one in the minimally coupled case. 
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1 Introduction 



The most attractive feature of zero temperature black holes is that they are the 
natural candidates as end-state of the evaporation process. Indeed, they represent 
the ideal setting where one can address the various issues connected to the quantum 
evolution of the black holes, such as for instance the problem of information loss 
(see e.g. gj). 

In spherically symmetric Einstein-Maxwell theory the only solution with this 
property is the extremal Reissner-Nordstrom (RN) black hole. Turning to the semi- 
classical theory, quantum corrections induced by the vacuum expectation value of 
the stress energy tensor due to matter fields modify the spacetime geometry and it 
is very important to check whether the resulting solution has still zero temperature 
or not. Perturbative corrections 0(h) to the classical geometry evaluated close to 
the horizon in four dimensions do not appear to answer unambiguously the above 
question 0, ||. It is clear that more information would come only if one knew the 
exact analytical solution to the semiclassical eqs. of motion. For the simple case of 
spherically reduced Einstein- Maxwell theory coupled with 2D minimal scalar fields 
Trivedi (3| was able to prove the existence of zero temperature solutions which re- 
duce, as h — > 0, to the extreme RN black hole. He also showed that although the 
energy density measured by an infalling observer close to the horizon diverges for 
the classical solution, the semiclassical configuration is regular there (only a mild 
singularity emerges in the second derivative of the scalar curvature ). The drawback 
of this analysis is that, due to the special type of matter fields used, these results do 
not have an obvious four dimensional interpretation. In order to improve this study, 
we consider here a more realistic 2D model that recently has received a lot of atten- 
tion . We employ a 2D conformal scalar field nonminimally coupled to the dilaton 
field, which classically corresponds to the s-wave sector of a 4D minimal scalar field 
(this model was first studied in 0). We will perform a numerical integration of the 
semiclassical equations of motion and show good evidence that zero temperature 
black holes exist in this theory. In particular , we inspect in detail the spacetime 
geometry in the region close to the horizon and compare with the results one gets 
by numerical integration of the minimally coupled case. 

The outline of this article is the following. In Section 2 we briefly review the 
spherically reduced Einstein-Maxwell theory and its zero-temperature solution, the 
extreme RN black hole. In Section 3 the matter model we shall use will be introduced 
and the expression of the (T a b) in the extreme RN background derived. In Section 
4 we numerically solve the backreaction equations and finally section 6 contains a 
discussion of our results and a comparison with the case analysed in || . 
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2 Einstein-Maxwell theory in D=2 

Let us start with Einstein-Maxwell theory in four dimensions 

S = Sq + Sem, (1) 
where Sq is the Einstein-Hilbert action [] 

S G = ±fd*xf^RW, (2) 

and Sem denotes the action associated to the electromagnetic field 

Sem = ~ I d'x^g^F 2 . (3) 

is the 4-dimensional scalar curvature and F 2 the strength of the electromagnetic 
field F^ v . Assuming spherical symmetry, the 4D metric can be written 

ds 2 = g [2) b dx a dx h + e -^ Xa) dn 2 , (4) 

(2) 

where g a l{x a ) (a,b — 1, 2) is the 2-dimensional metric in the (r — t)-plane, 4>(x a ) the 
dilaton and dfl 2 = d6 2 + sin 2 9d(p 2 the line element of the unit two-sphere. Dimen- 
sional reduction of the Einstein-Hilbert action (0) can be performed by integrating 
over the angles 9 and 4> 

S® = l -J d 2 x/^) e - 2 ^ (i2< 2 > + 2(V0) 2 + 2e 2 *) . (5) 

Proceeding similarly and considering F^ v = F^ v (x a ) the Maxwell action becomes 



o(2) 
°EM 



\ J d 2 x^a^e- 2 <f>F 2 , (6) 



where F 2 represents the field strength of a 2-dimensional gauge field. Black hole 
solutions of the theory defined by 

5 (2) = 4 2) + Sf M (7) 
are given by the Reissner-Nordstrom solution 

ds 2 = -f(r)dt 2 + J^dr 2 , e~^ = r 2 , (8) 

with 



3 We use units where h = G = c=ka = l 
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The parameter M is the ADM mass and Q the electric charge (the corresponding 
field strength is F rt = Q). The equation / = has two solutions for M > \Q\ given 
by r ± — M ± y/M 2 — Q 2 . r + = rh and r_ are, respectively, the event horizon and 
the inner horizon. The Hawking temperature is 

T„ = «!. (10) 
2nr 2 

Vanishing of Tjj, i.e. M = \Q\, defines the extremal configuration for which r + = 
r_ = Th = M. 



3 Matter fields 

In order to inquire on the existence of zero-temperature solutions in the semiclassical 
theory we must couple in eq. (0) to quantized free matter fields. In [f|] it was 
considered a 2D minimally coupled scalar field described by the classical action 



S M = --J d 2 x^J-g( 2 )(Vf) 2 (11) 
which, after quantization, yields the well-known Polyakov effective action 

This action can be formally obtained by functional integration of the trace anomaly 

(T) = —. (13) 

As it was mentioned in [Q, due to the 2d origin of this field the results one gets using 
Stat = + Sm + S e ff do not have an obvious four dimensional interpretation. To 
start with, we shall instead consider a 4D minimally coupled scalar field 



S { S = ~ / d'x^iVff. (14) 



16vr 

In a spherically symmetric spacetime, the matter fields can be expanded into spher- 
ical harmonics, the s-wave sector / of the scalar field / depending only on t and r. 
For the s-wave field / dimensional reduction gives the 2D action: 

S® = -~J *xy[^>e-»{Vf)*. (15) 



4 Usually, in order to make physical sense of the semiclassical approximation one considers N 
matter fields and consider the large N limit while keeping Nh fixed. In this way the quantum 
corrections due to the other fields can be neglected. 
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Comparison with the scalar field in ([□]) shows that the field /, though still 2D 
conformal, has acquired a nontrivial coupling with the dilaton field <fi. The corre- 
sponding trace anomaly has additional ^-dependent terms 



(T) 



247T 



R {2) -6(V0) 2 + 6D< 



(16) 



Performing a functional integration of this expression we get the following effective 
action §, 0, @ 



o(2) 



1 

24 



-g 



(2) 



1^(2)1^(2) _ I( V0 )2l i? (2) + I 0i? (2) 

48 □ 4 V ' □ 4 



(17) 



where the first nonlocal term is the same as in (|T2| ) . It is important to point out 
that unlike (|i~2"D this effective action is not exact. Unphysical results obtained for 
the evaporation of Schwarzschild black holes ||, suggest that, at least at finite 
temperature, Sfh must be modified by the addition of conformally invariant (local 



1C]. Considering instead zero-temperature configurations 



and nonlocal) terms | 
S^ff gives physically meaningful results [10J . 
In conformal gauge 

ds 2 = —fdudv (1£ 
this action becomes local (i.e.^R^ = —Inf) and for static configurations 



f(r) 



where 



u 



dr 

W) 



(19) 
(20) 



the components of the 2D stress energy tensor read 



/ T (2)\ 
\ uu J 



T 



(2) 



1 



96tt 



1 



64tt 



k 



ln/- 



k 



/7-(2) 
\ uv 



where the notation 



OZ7T 



96tt 



k 



V k" 1 

4+44' 



-2<f> 



+ 2- 

k 



k" 



(21) 
(22) 

(23) 

has been introduced and the prime denotes derivative with respect to r. Consid- 
ering now the dependence on the dilaton field, another relation can be deduced by 
functional differentiation of the effective action QT7D with respect to the dilaton 



This term is specific to the effective action considered and does not appear in the 
Polyakov theory. In a 4D viewpoint it is related to the tangential pressure (P) = 
T e e ) = (if) through the relation ([§], g) 

1 £9 (2) 

(P) = , (25) 



87re- 2 ^-c/( 2 ) S( f> 



For the particular case of the extremal Reissner-Nordstrom black hole / = (I — 
M/r) 2 , k = r 2 we obtain the following results (see also fTT|j ) 

where the first term on the r.h.s of these equations comes from the Polyakov con- 
tribution to the effective action. Also, from Eqs. (2J) and (|25| ) we obtain the 4D 
tangential pressure: 

^ ^ 167r 2 r 4 ( r)( r) ^ 167r 2 r 5 ( r ) ^ ^ 

Another important physical quantity is 

(gr-jjf) _ 4<t mm > 
/ / 2 ' 



(29) 



which is proportional to the energy density measured by an infalling observer |T3|| . 
Equation (EE) leads to 



1 M 1 

H n In 



M 
1 

r 



6?r r 2 (r-M) ' 2vrr 2 ~~ r ' ^ 

As in the minimally coupled case F diverges when r — > M. The term ~ 1/ (r — M) 
is the same as that found in the Polyakov theory 0], but despite its presence it 
was shown in M that the corresponding semiclassical zero-temperature solution is 
regular at the horizon (only a mild divergence is present in the second derivative 
of the scalar curvature R ). In our case , in addition to this term there appears a 
subleading logarithmic divergence ~ hx(r — M), which is present also in the analytic 
approximations in four dimensions proposed in JT2"| , as well as a nontrivial pressure 



(P) eq. (28) . As stressed in the first of Refs. , the divergence of F on the horizon 
of the classical extreme black hole makes the perturbative expansion in powers of H 
to break down there. The calculations performed in are motivated by the fact 
that F has been proven to be finite at r = M numerically in D=4 [13|, but due to 



the result (p0|) reliable near horizon calculations for zero temperature black holes 
performed using the effective action (|i~7l) must be nonperturbative in %. Similarly, 
the 0(h) results presented in |T4| (obtained by considering near-extreme black holes 
in the near horizon region) do not appear to have much physical meaning. 
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4 Backreaction 



We now come to the main question addressed in this paper: do self-consistent zero- 
temperature black holes exist in the semiclassical theory ? For this purpose, we 
need first of all to write down the semiclassical Einstein equations, which can be 
derived by differentiation of the action S^j = + Sff + S^X (see eqs. (0), (p~5|) 

_____ 

and fllTD ) with respect to the 2d metric g ab and the dilaton field 0. In conformal 
gauge (p~8|) and considering static configurations (|T9|) , (|20|) the relevant expressions 
concerning the matter part of the action have been derived in ([21]) , ( p2|) and (|24]). 
The corresponding quantities coming from the gravity and electromagnetic actions 
(where F 2 = —2^-) can be easily obtained by differentiation of (||) and (|6]). 
The uu (or vv ) constraint reads 



iiji 



U4) *+{(/*/-§</ 



t\2 



where the coefficient £ 



+3£ 



h 

12k 



(31) 



has been introduced (we have reintroduced H in the 



formulas in order to make the distinction between classical and quantum terms more 
clear) . The equation obtained by varying the trace of the metric (i.e. g uv ) reads 



= -2 + f'k' + fk" + 2^- + U" 

k 



3£ 



V k" 1 

4 + 4 - \f 



(32) 



Finally, differentiation with respect to gives 

! / // \ " 





J 2 \ k k k k 2 



f'k' f" . J f'k' 

-w-T +lnf [-w 



Ak'f f k' r 
k 3 k 2 



(33) 



For £ = (PH), (|32| ) and (^) are the classical eqs. of motion , for which the only 
zero-temperature configuration is the extremal Reissner-Nordstrom black hole / = 
(1 — M/r) 2 , k = r 2 . In the quantum terms, we have separated the ones multiplying 
£, coming to the Polyakov contribution to the effective action and present also in 
[01], and those proportional to 3£ representing the additional contributions in the 

(2) 

effective action S^/f (Eq. (|17D). 

As the three previous equations involve only two independent functions f(r) 
and k(r), one is redundant. Indeed, they are related through the Bianchi identities 
combined with the "nonconservation" eqs. for the matter part p| 



i(2)a 



we 



-20 



0. 



(34) 
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Also, as these non-linear differential equations involve the second order derivatives 
of / and k, two boundary conditions on these functions are required to determine 
them uniquely. For a zero temperature black hole, natural boundary conditions can 
be imposed on the function / at the horizon. First of all, / has to vanish there. 
Moreover, in the gauge used the temperature of the black hole takes the simple 
expression 

Tb = £, (35) 

where 

- = \ f\r=r h (36) 

is the surface gravity and denotes the radius of the horizon. So T# = means 
f = at r = r h . 

In the Polyakov case, starting from these boundary conditions Trivedi |4| has found 
the form of the exact solution (nonperturbative in H) close to the horizon as a (non- 
analytic) expansion in powers of the coordinate distance from the horizon r — rh- 
In our case, however, the terms proportional to In f complicate exceedingly this 
analysis and seem to prevent an expansion in closed form of the solution for small 
values of r — analogous to that proposed in || . 

A numerical resolution has then been undertaken and the boundary conditions 
have been imposed at infinity [15], where the solution is to a good approximation the 



extreme RN. In this region apart from a finite but very small renormalization of the 
classical mass (following || it is M R /Q = 1 + O(^)) the first quantum corrections 
to the spacetime metric are of the order 0(l/r 3 ). 

To start with, we have introduced dimensionless variables and functions in the dif- 
ferential equations that we have to integrate: 

x = r/Q, ~k(x) = k(r)/Q 2 , f{x) = f(r). 

This means choosing the black hole charge to be the natural unit of length. 

Our numerical integrations have been performed using the im-constraint equation 

H) and the 0-equation (|3~3"D for the value — = 1CT 5 . All along our calculations, 



Q 2 

the solutions of these equations have been checked to be compatible with eq. 
as well with a precision less than 10~ 7 . 

In order to probe the accuracy of our numerical simulations, we have first considered 
the integration of the semiclassical eqs. for the minimally coupled case (i.e. discard- 
ing the terms proportional to 3£ in Equations (|3lD , fl32|) and fl33|) ) and compared 
the numerical results with the form of the exact solution close to the horizon given 
by Trivedi We find that, for a value of the horizon xp — 1, 0229 (in units where 
Q = 1), i.e. with a deviation of about 2,3% from the classical value, the functions 
/ and /' behave like those in [|J with a precision of about 5 • 10 -4 when x — > xp and 
that k is accurate with a precision 5 • 10~ 5 and k' with the accuracy 2 • 10~ 4 . We can 
then expect the global precision of our simulation to be at least about 2 ■ 10 _2 %. 
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Considering now the nonminimally coupled case, we have integrated Equations ( ]5T| ) 
and (|33|). The results of this simulation are illustrated by the plots on the left of 
Figs.[l]-f| where the functions / and k and their first derivatives have been shown for 
x varying from the horizon xp> = I, 0378 to 5 (in units where Q — 1). To facilitate 
the comparison, the same functions in the minimal case (the plots on the right of 
Figs. [I|-^p have been reported for x varying from xp = 1, 0229 to 5. 

5 Discussion and conclusions 

Our numerical simulations presented in Figs. 1-4 and the comparison with the 
corresponding solution of the Polyakov theory appear to give good evidence that 
zero-temperature configurations exist in this theory. To get some insights from 
these results, we have compared the numerical values close to the horizon of these 
solutions with those obtained in the minimally coupled case. It turns out that the 
differences between the values of the functions / and /' for the two theories are less 
than 8 • 10 -4 (as deduced previously, the numerical precision is about 2 • 10~ 4 ). The 
same reasoning applies to the function k with a difference less than 8 • 10~ 5 . The first 
noticeable difference between the two theories appears at the first order derivative of 
the function k: the value of k! on the horizon for our model is estimated at —18, 094 
compared to +2, 075 for the Polyakov case. Going further in the derivatives of the 
function / we have that the third derivative of / blows up at the horizon for the 
nonminimally coupled field compared to the divergence of the fourth derivative in 
the minimal case. Moreover, it is interesting to stress that the value of the coordinate 
x at the horizon xp> = 1,0378 (in units where Q = 1) differs from xp = 1,0229 of 
about 1, 5% (up to a numerical precision of about 2 • 10 _2 %). 

Curvature invariants can be easily calculated starting from the results presented 
here. The 2D Ricci scalar = —f" is finite on the horizon , as it is shown in 
the plot of Fig. 5. Moreover, finiteness of k and kl are enough to prove that also 
the corresponding four dimensional scalar curvature is regular as the horizon 
is approached. A similar conclusion, despite the divergence of F at the horizon of 
the extreme RN black hole, was found in the minimally coupled case ||]. In our 
case, since the leading term in eq. ( |30| ) is the same as in the Polyakov theory, it 
is reasonable that the main conclusion about the regularity of the geometry at the 
horizon is unchanged. The difference with respect to the case analysed in || is 
that the "mild" singularity appearing at the horizon in the second derivative of the 
curvature close to the horizon (~ /"") is replaced by a "stronger" divergence in its 
first derivative (i.e. /"') . This is due to the logarithmic divergence ln(/) appearing 
in eq. (|30|) as well as in eqs. ([31]) and (|33f) . 

In conclusion, by numerical integration of the two-dimensional semiclassical 
equations of motion for the case of spherically reduced Einstein-Maxwell theory 
and a scalar field nonminimally coupled to the dilaton we have found solutions de- 
scribing zero-temperature black holes. Similarities and differences with respect to 
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the simpler minimally coupled case have been studied in detail. Due to the intrinsic 
four dimensional nature of the matter field used, our results could be relevant in or- 
der to address the same issue in the physical world D = 4 . Finally, following Refs. 
[0] an interesting extension of this work would be to check whether the solution 
found here does indeed represent the end-point of the evaporation process (for the 
minimally coupled case this problem has been addressed both in the near-horizon 
approximation [[U| and in the whole spacetime numerically [|17]]) . 



Acknowledgements 

We would like to thank R. Balbinot and P. Nicolini for useful discussions and for 
collaboration at an earlier stage of this work. C.B. was supported by the Italian 
Minister of Foreign Affairs. 



10 




Figure 5: Values of /" = 



— for our model. 
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